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(54) Soft value calculation for multilevel signals 

(57) A sub-optimal method is disclosed for calculat- 
ing the reliability values (soft values) for the bits of a mul- 
tilevel signal. The log-likelihood values are approximat- 
ed using only the dominant terms, so called max -log ap- 
proximation, that is for each bit position only the two 
closest signal symbols of opposite bit value (S 8 ,S 6 ) are 
considered in the sum. The used modulation scheme is 
16-QAM together with Gray-labelling. Two versions of 
approximation are proposed: one version consists of us- 
ing the two distances between the received value and 



the two closest symbols of opposite bit value (S^o^). In 
order to simplify and speed up the calculation, the sec- 
ond version consists of using the distance between the 
two closest symbols (5 3 ) to approximate the distance 
between the second closest symbol and the received 
value. Furthermore, precalculated results are stored in 
look-up tables to speed up the calculation. 

Possible applications are especially bit interleaved 
coded modulation (BICM) together with soft-input de- 
coding. It is also of interest for TCM and BCM schemes. 



0111 


0101 


1101 


1111 


o 

s 1 


o 

s 2 


o 

s 3 


o 

s 4 


0110 


0100 4 


5 3 1100 


1110 


o 

S 5 




o 










x — 201 


0010 


0000 


1000 


1010 


o 

s fl 


o 

s 10 


o 

8„ 


o 

s 12 


0011 


0001 


1001 


1011 


o 

Si3 


o 

8l4 


o 

8« 


o 



Fig. 2 



Printed by Jouve. 75001 PARIS (FR) 



1343285A1 I > 



EP 1 343 285 A1 



Description 

[00011 This invention relates to digital communications systems and, more particularly, the generation of soft reliability 
values for multilevel signals. . 

5 [0002] Within the field of digital communications, multilevel modulation is used to map a number of bit sequences to 
a signal alphabet comprising a number of signal symbols, i.e. a number of points in signal space. For example, a bit 
sequence may be mapped onto a point in a complex signal space. A signal alphabet of size M allows log 2 (M) bits to 
be mapped to each symbol. However, when symbols are received at a receiver, they may be affected by noise, thereby 
affecting the decoding of the signal when retrieving the transmitted bit sequence. If multilevel modulation is used in 

10 conjunction with channel coding, many channel decoders, such as iterative decoders based on the BCJR algorithm, 
require so-called soft bit values as an input. A soft bit value corresponds to a reliability value of a single bit being 0 or 1 . 
[0003] Examples of multilevel modulation include multi-amplitude level modulation in Pulse Amplitude Modulation 
(RAM), multi phase level modulation in Phase Shift Keying (PSK), multi signal point modulation in Quadrature Amplitude 

Modulation (QAM). , (f . 

15 [0004] For example, an emerging technology for wideband digital radio communications of Internet, multimedia, 
video and other capacity-demanding applications in connection with the third generation of mobile telephone systems 
is the evolving Wideband Code Division Multiple Access (WCDMA) specified as part of the 3GPP standardisation 
organisation. Within this technology, High Speed Downlink Packet Access (HSDPA) Is provided including a high speed 
downlink shared channel (HS-DSCH) which uses 16QAM. In 16QAM for example, M=16. i.e. each symbol in the signal 

20 alphabet represents 4 bits. Future releases may comprise even larger constellation sizes such as 64 QAM. 

[0005] It is known how to convert signal symbols to soft bit values by calculating all distances in signal space between 
the received symbol and all signal points of the signal alphabet. In particular, in order to obtain optimal performance, 
a likelihood ratio is calculated depending on corresponding jsums of probabilities where the probabilities are functions 
of the calculated distances. It is further known that in the calculation of a likelihood ratio the sums of probabilities may 

25 be approximated by the dominant contributions to the sums of probabilities in the likelihood ratio (A.J. Viterbi, "An 
intuitive justification and a simplified implementation of the MAP decoder for conventional codes", IEEE Journal on 
selected areas in communications, 16(2), February 19g8). 

[0006] Even though this approximation significantly reduces the computational complexity while only causing a neg- 
ligible loss in performance, it is a problem of the above prior art method that it still requires a calculation of all distances 

30 to all the signal points in order to determine which two are actually needed for the calculation of the likelihood ratio. 
For example, in a 1 6QAM modulation, 1 6 distances have to be calculated for each received symbol. In particular, if a 
high rate of symbols need to be decoded, e.g. several hundred symbols per millisecond, the above problem is partic- 
ularly severe. ^ 4 
[0007] The above and other problems are solved when a method of generating a reliability value for a received 

35 multilevel signal in relation to a number of predetermined signal symbols each associated with a corresponding bit 
sequence including a first bit position; the reliability value being indicative of likelihood information of receiving said 
multilevel signal is characterized in that the method comprises the steps of 

- identifying a first one of the number of signal symbols as being closest to the received multilevel signal; and 
40 - estimating the reliability value based on stored information related to a second one of the number of signal symbols 
being closest to the first signal symbol and corresponding to a different binary value at the first bit position of the 
respective associated bit sequence than the first signal symbol. 

[0008] Consequently, according to the invention, the closest signal point needed for calculating the likelihood ratio 
45 is determined first, such that only the corresponding distances to the identified signal points need to be determined. 
In this way, the computational complexity is reduced significantly, since not more than two distances need to be cal- 
culated in order to determine the likelihood ratio for each bit. Hence, the above method is particularly well-suited for 
low-complexity implementations in mobile receivers, as it reduces the required computational resources. 
[0009] According to a preferred embodiment of the invention, the stored information comprises an identification of 
50 the second signal symbol, and the step of estimating the reliability value further comprises the steps of 

- determining a first distance between the received signal and the first signal symbol; and 

- determining a second distance between the received signal and the second signal symbol. 

55 [0010] Hence, given a signal constellation, as for each signal symbol and for each bit position it is known which one 
of the other signal symbols having an opposite value at that bit position is closest to that signal symbol, this information 
may be stored. Therefore, by a simple look-up operation the closest signal symbol may be identified and the distance 
between the received signal and the closest signal symbol may be calculated. 
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[0011] According to another preferred embodiment of the invention, the stored information comprises a second dis- 
tance between the first signal symbol and the second signal symbol, and the step of estimating the reliability value 
further comprises the step of determining a first distance between the received signal and the first signal symbol. 
Hence, alternatively or additionally to storing an identification of the signal symbol with opposite bit value at a given 

5 bit position being closest to the first signal symbol, the actual distance between the first and the second symbol may 
be pre-calculated and stored. Once the first signal symbol is identified, this distance may be looked up and used as 
an approximation for the distance between the received signal and the second signal symbol. Consequently, a further 
reduction in computational complexity is achieved, since only one distance has to be calculated. 
[0012] When the step of estimating the reliability value comprises the step of determining a polynomial function of 

10 the first distance and the second distance, multiplied by a predetermined constant, the computational complexity is 
further reduced, as no logarithm needs to be calculated. In one embodiment, the polynomial function is a difference 
of the squared distances. 

[0013] According to yet another preferred embodiment of the invention, the stored information is indicative of a 
number of functional relations between the received multilevel signal and the reliability value, and the step of estimating 

15 the reliability value further comprises the step of selecting a functional relation of said number of functional relations 
dependant on the first signal symbol and the first bit position. Hence, the calculation of the likelihood value only com- 
prises the step of calculating the corresponding stored function. Preferably, the functional relationship isalinearfunction 
of a signal component, thereby reducing the calculation to a multiplication operation and an adding operation. 
[0014] According to yet another preferred embodiment, the stored information comprises, for each signal symbol 

20 and bit position an approximation of the corresponding reliability value. Hence, according to this embodiment, an ap- 
proximation of the reliability value may be directly looked up once the closest signal symbol is identified, thereby pro- 
viding a computationally very, efficient method which eliminates the need of online distance calculations. 
[0015] Preferably, the stored information is stored in a look-up table indexed by the number of signal symbols and 
the bit positions of the number of bit sequences, thereby providing fast access to the information. 

25 [0016] In one embodiment of the invention, the method further comprises the step of providing the reliability value 
as an input to a decoder, e.g. an iterative decoder using the BCJR algorithm or any other decoder using soft values 
as an input. It is an advantage of the invention that it provides an accurate and resource-efficient approximation of soft 
values as an input to such decoders. 

[0017] The first signal symbol may be identified by comparing the signal components with predetermined thresholds 
30 or decision boundaries, for instance by means of a slicer, i.e. a circuit which compares a signal with predetermined 

thresholds. Hence, a fast and computationally inexpensive method is provided for identifying the closest signal symbol 

without the necessity of calculating all distances between the received value and all signal symbols. 

[0018] It is a further advantage of the invention that cost-effective, standard components may be employed when 

implementing a method according to the invention. 
35 [0019] When the likelihood information comprises a log-likelihood ratio, a high performance quality is achieved, as 

the use of a log-likelihood corresponds a theoretically optimal way of calculating soft reliability values. However, other 

methods of calculating likelihood information may be employed, such as a log-likelihood of the signal power. 

[0020] In a preferred embodiment of the invention the step of identifying the first signal symbol as being closest to 

the received multilevel signal comprises the step of identifying the first signal symbol as being closest to the received 
40 multilevel signal with respect to a Euclidean distance measure in a signal space, as the Euclidean distance is directly 

related to the probabilities of a likelihood calculation. Alternatively, other suitable known metrics may be used instead 

Euclidean distances. 

[0021] The signal space may be a real or complex signal space. For example, in QAM modulation two amplitude- 
modulated signals are transmitted on a single carrier, but shifted in phase by 90 degrees. Hence, the resulting signal 
4* points may be represented in the complex plane representing the so-called in-phase (I) and quadrature (Q) components 
of the QAM signal. 

[0022] When the number of signal symbols are associated with the number of bit sequences such that the bit se- 
quences associated with all nearest neighbours of each signal symbol only differ from the bit sequence of that signal 
symbol at one bit position, the error rate of the transmission system is reduced. This form of mapping is referred to as 
so Gray mapping. 

[0023] . The invention further relates to an arrangement for generating a reliability value for a received multilevel signal 
in relation to a number of predetermined signal symbols each associated with a corresponding bit sequence including 
a first bit position; the reliability value being indicative of likelihood information of receiving said multilevel signal; char- 
acterised in that the arrangement comprises 

55 

- first processing means adapted to identify a first one of the number of signal symbols as being closest to the 
received multilevel signal; 

- storage means adapted to store information related to the first signal symbol and a second one of the number 
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of signal symbols being closest to the first signal symbol and corresponding to a different binary value at the first 

bit position of the respective associated bit sequence than the first signal symbol; and 

- JSi ^processing means adapted to estimate the reliability value on the basis of the stored information. 

5 [0024] The arrangement may be implemented by any processing unit, e.g. a programmable micro P™f^ an 
a P p¥cition-specificintegratedcircuit,oranotherintegratedc^ 

comprises general- or special-purpose programmable microprocessors, Digital Signal Process'* > (DSP ) Applet on 
Specif c integrated Circuits (ASIC), Programmable Logic Arrays (PLA). Field Programmable Gate ^(<^M*j 
or a combination thereof. The processing means may be a CPU of a computer, a m.croprocessor, a smart card, a SIM 
to card or the like. The first and second processing means may be separate processing means, ejj. separate circuits 

[oOMr^term storage means includes magnetic tape, optica, disc, digital video disk (DVD), compact dtec (CD or 
CD-ROM), mini-disc, ha* disk, floppy disk, ferro-electric memory, electrically e-sab.e programmab J-ad or^y mem 
orv (EEPROM) flash memory. EPROM. read only memory (ROM), static random access memory (SRAM), dynamic 
Z^ ^^^iO^M), synchronous dynamic random access memory (SDRAM), ferromagnetic memory, 
optical storage, charge coupled devices, smart cards, etc. 

[0026] Furthermore, the above discussed features and steps of the method according to the invention may be incor- 
porated in the above arrangement according to the invention. 

20 10027] The invention further relates to a device for receiving multilevel signals compris.ng an arrangement as de- 
scribed above and in the following. , . . „, 
[0028] The device may be any electronic equipment or part of such electronic equipment, where the term electronic 
KLnt incites computers, such as stationary end portable PCs. stationary and portable ^ 
equipment. The term portable radio communications equipment includes mobile radio torm.nals such as mob.le tele- 

ss nhonas Daaers communicators, e.g. electronic organisers, smart phones, PDAs, or the like. 

[0020?' SentrwS explained more fully below in connection with preferred embodiments and with reference 
to the drawings, in which: 

figs. 1a-b schematically show a receiver according to an embodiment of the invention; 

30 

fig. 2 shows an example of a signal constellation with 16 signal symbols; 
fig. 3 shows a flow diagram of a method according to a first embodiment of the invention; 
35 fig. 4 shows an example of a look-up table according to the embodiment of fig. 3; 

fig. 5 shows a flow diagram of a method according to a second embodiment of the invention; 
fig. 6 shows an example of a look-up table according to the embodiment of fig. 5; 

40 

fig. 7 shows another example of a signal constellation with 16 signal symbols; 

fig. 8 shows a flow diagram of a method according to a third embodiment of the invention; and 

49 fig. 9 shows an example of a look-up table according to the embodiment of fig. 8. 

[0030] Fig. 1 a schematically shows a receiver according to an embodiment of the invention receiving a radio signal 
from a transmitter via a communications channel. Transmitter 101 is adapted to send a signal s via a noisy channel 
102 to receiver 103. The signal s represents one of a set of M signal points Si ... S M in a signal space where each 

so signal point is related to a respective bit sequence of log 2 (M) bits. In the presence of noise in the transmission , channe 
102 the receiver 103 receives a signal f which deviates from the transmitted signal s. In one embodiment the signal 
is a Code Division Multiple Access (CDMA) signal using a spread spectrum technique. The receiver 1 03 compnses a 
receiver circuit 1 07 for transforming the received spread spectrum signal into the signal symbol r. The receiver further 
comprises a channel decoder 106 for decoding the received signal symbol r, e.g. a BCJR or Viterbi decoder. The 

55 decoder 106 requires soft bit values as an input. Hence, the receiver 103 further comprises a circuit 104 which is 
adapted to calculate soft values for the log 2 (M) bits of the received signal symbol r and to provide the calculated soft 
values to the decoder 106. According to the invention, the receiver 1 03 further comprises a memory 1 05, such as on- 
chip mem07, EPROM, flash memory, or the like, in which a look-up table is stored for use in an efficient calcuiat.cn of 
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the soft values by the circuit 104, preferably as described in connection with figs. 3-9. 

[0031] Fig. 1b schematically shows a more detailed block diagram of the receiver circuit 1 03 of fig. The circuit 107 
comprises a RAKE receiver 110 suitable for receiving CDMA signals, i.e. a receiver which uses several baseband 
correlators to individually process several signal multipath components. The correlator outputs are combined to achieve 

5 improved communications reliability and performance (see e.g. "Digital Communications" 4th Edition, by John G. 
Proakis, McGraw-Hill, 2000). The sampled received radio signal r* is fed to the RAKE receiver 110 which generates 
the signal symbol r to be decoded. The circuit 107 further comprises a channel estimator 111 and a noise estimator 
112, e.g. implementing any suitable channel estimation and noise estimation technique known in the art. The channel 
estimator receives the received radio signal r\ identifies up to N different radio paths or channel taps and estimates 

10 corresponding delays A k , k=1 ,...,N, andcomplex channel estimate h r = (h r1 ... ,h rN ) of these paths. The channel estimator 
111 further provides a set of complex combiner weights w = [W n W 2 ... W N ] T to be used by the rake receiver. Here, [] T 
denotes a transposed vector. For example, the weights may be determined according to an optimisation criterion, such 
as maximising the received signal energy. 

[0032] The calculated delays A k and the combiner weights are provided to the rake receiver 1 1 0. The RAKE receiver 
15 no comprises delay circuits 1 1 5 which delay the incoming signal according to the N channel taps. Further, the receiver 
110 comprises circuitry 116 for multiplying the N delayed versions of the received signal with a spreading code c for 
dispreading the spread spectrum signals and circuitry 11 7 for summing the signals to fonm a radio symbol. Furthermore, 
the rake receiver 110 comprises multiplier circuitry 118 for multiplying each of the N radio symbols with the combiner 
weights (w k ) # , k=1 ,...,N, where ()* denotes complex conjugation. Finally, the RAKE receiver 110 comprises an adding 
20 circuit 119 which combines the weighted symbols to form the received symbol estimate r which is fed to the soft value 
calculation circuit 1 04. 

[0033] When using a multilevel constellation of signal points S 1 ...S M the amplitude information should be maintained 
in order to ensure successful demodulation in the receiver. Consequently, the reference points S 1 ... S M should be 
scaled properly. In the following it is assumed that the channel estimator 111 estimates the channel gain on the basis 
25 of a reference channel h r which has a channel gain that may be different from the actual gain of the traffic channel, e. 
g. a HS-DSCH. The gain difference between the reference channel and the traffic channel may be denoted with g. 
Hence, the received symbol r after the RAKE receiver 110 may be expressed as 

30 r=gw H h r s + n, 

where is the Hermitian conjugate of w, w H h r denotes an inner product, s is the transmitted symbol and n is a noise 
term, e.g. representing additive white Gaussian noise (AWGN). The gain parameter g is signalled to the receiver, w is 
selected by the combiner in the receiver, and h r are the channel estimates. Hence, at the receiver, the reference signal 
35 symbols S 1 ... S M may be scaled appropriately, according to 

S j = gw H h r S j( j=1 M. (0) 

40 [0034] In fig. 1 b, the receiver circuit 1 07 comprises circuit 1 1 3 adapted to calculate the above scaling factor g w^n, 
and a multiplier circuit 114 for multiplying the reference symbols with the scaling factor, resulting in properly scaled 
signal symbols S t S M , which are fed into the soft value calculation circuit 104. 

[0035] Finally, the noise estimator 1 1 2 provides an estimate of the signal noise level a which is fed into the soft value 
calculation circuit 104. 

45 [0036] According to the invention, the soft value calculation circuit locates the signal symbol which is closest to the 

received signal r and calculates corresponding soft values L m for bit m, m=1 log 2 (M), e.g. according to one of the 

embodiments discussed in connection with figures 3-9. 

[0037] It is noted that the receiver circuit described in connection with figs. 1 a-b merely serves as an example, and 
the scope of the invention is not limited to the type of receiver, nor to the above scaling of signal symbols. 

so [0038] Fig. 2 shows a signal constellation with 16 signal symbols. The signal constellation comprises M=16 signal 
points S 1 through S 16 in a two-dimensional signal space, e.g. the l/Q components in a 16QAM signal constellation. 
Preferably, the signal points are distributed regularly, such that the distance to the nearest neighbours of each signal 
point is the same. The reference points may take values that suit the implementation in question. However, alternatively, 
other signal constellations may be chosen. In fig. 2, 16 different bit sequences 0000 through 1111 , each consisting of 

55 log 2 (M)=4 bits, are mapped onto the signal points S r S 16 . Preferably, the mapping of the bit sequences to the signal 
points is chosen such that the bit sequence of each signal point only differs from those of the nearest neighbours by 
one bit, thereby optimising the decoding performance. For example, in fig. 2, signal point S 8 has three nearest neigh- 
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bours, S 4 , S 7 , and S 12 . The bit sequence of S 4> i.e. 1111 , differs from the sequence 1110 of S 8 only at bit position 4, 
etc. Alternatively, other mappings may be chosen. 

[0039] For every bit position m mapped on a signal point, the signal points in the constellation may be divided into 
two sets where the signal points in each set have bit value 0 and 1 , respectively, at that position. In the following, the 
5 set of signal points with a 0 at the m-th position is denoted Ao <m , and the corresponding set with a 1 at the m-th position 
is denoted A, m . For example, in the example of fig. 2 and for m=1 , A n t1 ={S 3 , S 4> S 7 , S 8 , S 1t , S 12 , S 15 , S 16 } and Aq <,= 
{S 1( S 2 , S 5 , Se, S 9 , S n0 , S 13 , S 14 ). The sets are of equal size with M/2 elements each. 

[0040] When one of the signals S 1t ...,S 16 is transmitted over a noisy channel, the received signal will differ from the 
transmitted signal according to a corresponding distribution. The actual shape and width of the distribution of received 
10 signals depends on the characteristics of the noise. In fig. 2, the cross 201 represents an example of a received signal r. 
[0041] Prior to providing the received 1 6QAM radio symbols to a decoder, e.g. a turbo decoder, they are converted 
into soft values. Hence, a soft value is calculated for each bit of every 1 6QAM symbol. A soft value of the m-th bit in 
the sequence mapped to r may be defined as 
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, , P(s m =1lr) P(s m =1)P(r ls m =D ,_ P(rls m =1) 

m ' 109 P(s m =0lrf ,O9 P(s m =0) P(rls m =0f ,og P(r ls m =0) { } 



where s m is the m-th bit in the bit sequence represented by the transmitted signal, and P(s m =il r), 1=0,1, are the a 
20 posteriori probabilities of the bit s m where r is the received signal. It is noted that the second equality assumes that 
s m =1 and s m =0 are equally probable in the chosen alphabet. Otherwise, the overall ratio of probabilities should be 
taken into consideration in the following. However, this would only give rise to a constant factor. Hence, L m corresponds 
to a log-likelihood ratio of probabilities. The probabilities P(r ls m =i) in eqn. (1) may be written as 



P ( r |Sm=0=^ Z p (r.s). i = 0,1. (2) 



SeA it _ 



[0042] Hence, the calculation of the above probability involves a summation over W2 terms each including a joint 
probability P(r,s). This is a computationally expensive task, especially if M is large, e.g. M=64. 
[0043] In many applications, the above soft values may be approximated by 

max P(r,s) (3) 
Lm=, ° 9 max P(r,sr l ° g W^)' 

SGA 0,m 



where s, m i=0,1, are the signal points that result in the largest contribution to the sums in eqn. (2). Hence, in the 
calculation of the probabilities, the sums over M/2 terms are approximated by the their respective dominant terms, 
45 according to 

log £P(r,s)«log max P(r,s) (4) 

50 SeA », m S€A i.m 

[0044] The above approximation is often referred to as the "max log MAP" approximation which yields a good ap- 
proximation in cases where the above sums are dominated by one term, as for example in the case of Gaussian noise 
55 when the signal to noise ratio (SNR) is large. The above probabilities depend on the distances between the received 
signal r and the respective signal points. For example, in the case of additive zero-mean Gaussian noise with variance 
o 2 , the log-likelihood ratio of eqn. (3) may be expanded as 
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a 2 exp(-|r-S 1m | 2 /a 2 J 

Lm=log— f- ^"^ = a (l r - S 0.m| -|r-Hm| } (5) 

o - 2 exp[-|r-S 0)m p/a 2 J v J 

[0045] It is noted that it is assumed that the Sj m are scaled corresponding to the received signal according to equation 
10 (0) above. In the following, we define dj m =lr-Sj m l, for i=0,1 , to be the distances between the received signal r and the 
closest signal points in the sets \ m , respectively. For example, in fig. 2, for m=1 , d 0 «, corresponds to the distance 5 2 
between r and the closest signal point in Aq 1t i.e. S 6 , while d 1 A corresponds to the distance 5 1 between r and the 
closest signal point in A 1 A , i.e. S 8 . According to the invention, the likelihood ratio in equation (5) is obtained by first 
identifying the closest signal point S 8 , and then determining the distances 5-, and 6 2 , as will be described in greater 
is detail in connection with figs. 3-4. Hence, according to the invention, a computationally expensive calculation of all the 
distances between r and all the signal points S V ..S 16 in order to identify the shortest distances 5 n and is avoided. 
Alternatively to Identifying S 6 by means of a look-up table and then calculating 5 2 , the distance 63 between S 6 and S 8 
may be looked up and used as an approximation instead of 6 2> thereby saving additional computational resources. 
This will be described in greater detail in connection with figs. 5-6. A further embodiment of the invention will be de- 
20 scribed in connection with figs. 8-9. 

[0046] It is noted that, preferably, in the above estimation of the reliability values, a proper scaling of the signal points 
in the QAM constellation is taken into consideration. If this scaling is taken into consideration, the above log-likelihood 
ratio may be written as 
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Lm=a- 2 {|r- So ,m| 2 -|r-Hm| 2 } < 6) 



[0047] As described in connection with fig. 1b, when using a multilevel constellation as in the example of fig. 2, the 
amplitude information should be maintained in order to ensure successful demodulation in the receiver. Consequently, 
the reference points S jt j=1 ,...,16 should be scaled properly. If this scaling is taken into consideration, the above log- 
likelihood ratio may be written as 

L m = K.[|r-s 0im | 2 -|r-? 1>m | 2 ], (7) 

[0046] i.e. with the properly scaled signals 

r 



r = 



gw H h r 



3 it m=-%- m = 1 log2(M),t=0,1, (8) 

gw n h r 



and where 



K=(gw M h r ) 2 /o 2 (9) 



is a constant which depends on the signal to noise ratio. 
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bours. S 4> Sy, and S 12 . The bit sequence of S 4 , i.e. 1111 . differs from the sequence 1110 of S e only at bit position 4. 
etc. Alternatively, other mappings may be chosen. 

[0039] For every bit position m mapped on a signal point, the signal points in the constellauon may be divided into 
Jwo sets where the signal points in each set have bit value 0 and 1. respectively, at that Position. In the flowing he 
5 set of signal points with a 0 at the m-th position is denoted Ao. m . and the corresponding set with a 1 at the m-th posihon 
is denoted A, m . For example, in the example of fig. 2 and for m=1 , A, ,={S 3 . S 4 , S 7 . S 8 , S t „ S 12 , s 15 , & 16 ) ana 
/S. S« S= sL So Sin S11. S,a). The sets are of equal size with M/2 elements each. 

00401 When one of the signals S, S 16 is transmitted over a noisy channel, the received signal will differ from the 

transmitted signal according to a corresponding distribution. The actual shape and width of the distnbut.cn of received 
10 signals depends on the characteristics of the noise. In fig. 2. the cross 201 represents an example of a received s.gnal r. 
[00411 Prior to providing the received 16QAM radio symbols to a decoder, e.g. a turbo decoder, they are converted 
into soft values. Hence, a soft value is calculated for each bit of every 16QAM symbol. A soft value of the m-th bit .n 
the sequence mapped to r may be defined as 
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P(s m =1lr) P(s m =1)P(rls m =1) P(r ls m =1) 
L m = '°9 P(Sni=0 lr) = log P(s m =0)P(rls m =0) =l09 P(r ls m =0) 



0) 



where s m is the m-th bit in the bit sequence represented by the transmitted signal, and P(s m =.l r), .=0,1. are the a 
20 posteriori probabilities of the bit s m where r is the received signal. It is noted that the second equa^y assumes that 
s =1 and s m =0 are equally probable in the chosen alphabet. Otherwise, the overall ratio of probab.l.t.es should be 
taken into consideration in the following. However, this would only give rise to a constant factor. Hence, L m corresponds 
to a log-likelihood ratio of probabilities. The probabilities P(r ls„,=i) in eqn. (1) may be wrttten as 



P(r|s m =0=^ ZKr.s). i = 0,1. (2) 



seA u« 



[0042] Hence, the calculation of the above probability involves a summation over M/2 terms each includmg a joint 
probability P(r.s). This is a computationally expensive task, especially if M is large, e.g. M=64. 
[0043] In many applications, the above soft values L„ may be approximated by 

max P(r,s) (3) 

seA 1m . P(r|S 1m ) 

SeA 0.m 

where s, 1=0,1, are the signal points that result in the largest contribution to the sums in eqn. (2). Hence, in the 
calculation of the probabilities, the sums over M/2 terms are approximated by the their respective dominant terms, 
45 according to 



log £ p ( r ' s )* log max P(r,s) (4) 



so $€A j,m 5€A i,rn 



t0044] The above approximation is often referred to as the "max log MAP" approximation which yields a good ap- 
proximation in cases where the above sums are dominated by one term, as for example in the case of Gaussian noise 
55 when the signal to noise ratio (SNR) is large. The above probabilities depend on the distances between the received 
signal r and the respective signal points. For example, in the case of additive zero-mean Gaussian noise with variance 
o* the log-Nkelihood ratio of eqn. (3) may be expanded as 
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bits. Hence, it is a further advantage of this embodiment that it requires little storage capacity. In an embodiment where 
the resolution of the pre-computed distances is higher than log 2 (M) bits, i.e. ii>log 2 (M), processing time is traded for 
memory space in comparison with the embodiment of figs. 3-4. 

[0056] It is noted that additional storage space may be saved by only storing each distance once, i.e. in case the 
5 same distance appears in two or more entries of the table, a reference to that distance may be stored in one of the 
entries, instead. 

[0057] Alternatively, other layouts of a look-up table may be used. For example, in one embodiment, the look-up 
table may comprise all M(M-1 )/2 mutual distances between the signal points S k> thus requiring tiM(M-1 )/2 bits of stor- 
age. However, for M>4 this embodiment requires larger storage capacity than the embodiment of fig. 6. 
io [0058] Fig. 7 shows another example of a signal constellation with 1 6 signal symbols. As in fig. 2, the signal constel- 
lation comprises M=1 6 signal points S 1 through S 16 in a two-dimensional signal space, e.g. the l/Q components in a 
1 6QAM signal constellation. The signal points are distributed regularly, such that the distance to the nearest neighbours 
of each signal point is the same. In the example of fig. 7, they are assumed to be selected such that 

15 S k = x k +j-y k , where x k , y k G [-3d, -d, d, 3d], k=1 ,...,M, 

where d is an arbitrary constant and where j 2 = -1 . For example, d may be chosen to d=1 . However, alternatively, other 
signal constellations may be chosen. 
20 [0059] In fig. 7, 16 different bit sequences 0000 through 1111 , each consisting of log 2 (16)=4 bits, are mapped onto 
the signal points S 1 S 16 . Preferably, the mapping of the bit sequences to the signal points is chosen to be a Gray 
mapping, i.e. such that the bit sequence of each signal point only differs from those of the nearest neighbours by one 
bit, thereby optimising the decoding performance. 

[0060] As above, for every bit position m mapped on a signal point, the signal points in the constellation may be 
25 divided into two sets Aq m and A 1 m , where the signal points in each set have bit value 0 and 1 , respectively, at that 
position. 

[0061] Fig. 8 shows a flow diagram of a method according to a third embodiment of the invention. Again, this em- 
bodiment utilises the approximation of equation (7) for the calculation of the soft values L m . As in the embodiment of 
fig. 3, in the initial step 801 , a signal r is received. The received signal may be written as r = Re(r) + j lm(r) and, in the 

30 following the magnitude of the I- and Q components of r will be denoted by a = IRe(r)l and b = llm(r)l, respectively. 
After the received symbol is combined in the combiner, in step 802, the signal point s from the set of signal points S v .. 
S M , which is closest to r is identified. In this embodiment it is assumed that the constellation of signal points corresponds 
to the constellation of fig. 7. In fig. 7, each signal point corresponds to a decision region where the decision regions 
are separated by a set of decision boundaries 701 through 706. Hence, the closest signal point $ may be found by 

35 performing two comparisons of the inphase component and the quadrature component, respectively. For example, if 
Re(r) < 0 (decision boundary 705) and Re(r) <-2d (decision boundary 706) and if im(r)>0 (decision boundary 702) and 
lm(r)>2d (decision boundary 701), the received signal lies in the decision region corresponding to S v i.e. $= S n is the 
closest signal point. Subsequently, for each bit positions m= 1 , . . . , log 2 S (1 6)=1 .... ,4, the soft value L m may be calculated 
using the approximation of eqn. (7), assuming proper scaling. Consequently, in this example the soft value L 1 for the 

40 first bit is 

L t (S,) = K(l-a+j-b - (d+3j-d)1 2 - l-a+j-b - (-3d+3j-d)l 2 ) = K (8ad-8d 2 ), 

45 [0062] Hence, in the above equation, instead of computing two distances squared, each involving a calculation of 
the type lx+jyl 2 , the soft value may be calculated by scaling the inphase amplitude a of the received symbol with 8dK 
and, subsequently, by adding a constant -SKd 2 . It is further noted that the constant d may be chosen as any suitable 
positive real number. 

[0063] The remaining three soft values for a received signal in the decision region corresponding to S 1 are accordingly 

50 

L 2 (S 1 ) = K(l-a+j-b - (-3d+3j.d)l 2 - l-a+jb - (-3d-j-d)l 2 ) = K (-8bd+8d 2 ) 
55 L Z (SJ = K(l-a+j-b - (-d+3j-d)l 2 - l-a+j-b - (-3d+3jd)l 2 ) = K (4ad-8d 2 ) 

M s i) = KO-a+j-b - (-3d+3j-d)l 2 - l-a+j b - (-3d+jd)l 2 ) = K (4bd-8d 2 ), 
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as in the constellation of fig. 7 the closest symbols with opposite second, third, and fourth bit compared to are S 9 
= -3d-jd. So = -d + 3jd, and S 5 = -3d + jd, respectively. 

[0064] The table 808 of fig. 9 illustrates the calculated soft values for all decision regions corresponding to the symbols 
Si S 16 and for all bits. m=1 4. As can be seen from table 808, all soft values may be calculated by scaling one 

5 of the inphase component a or quadtrature component b of the received signal r and subsequently adding a constant. 
Hence using the pre-calculated equations of table 808, the soft values may be calculated in a very efficient way. In 
one implementation, each entry of the look-up table 808 may comprise the scaling factor, the constant to be added 
and a bit indicating whether it is the inphase component a or the quadtrature component b of the received signal r 
which is to be scaled for a given soft value. Preferably, the table is indexed by the decision region and the bit values. 

10 It is noted, however, that many of the entries of table 808 are identical. Consequently, it will be apparent to a skilled 
person that table 808 may be stored in a memory efficient manner, e.g. by storing a list of the distinct entries and, in 
table 808 referring to the corresponding list members. In general, it is noted that constellations which are Gray coded 
or show another regularity, the redundancy of the entries in table 808 may be utilised to reduce the memory consumption 
of table 808. 

15 [0065] Referring again to fig. 8, in steps 804-805 the soft values for the identified decision region and for all bits are 
calculated In step 804, the relation to be calculated, i.e. the scaling factor and the constant to be added, are retrieved 
from a stored table 808 In memory, e.g. a look-up table as shown in fig. 9. The retrieved relation is calculated in step 
805 resulting in the soft value for the corresponding bit number 

[0066] It is noted that the processing load in the receiver may further be decreased by pre-calcuiating the relations 
20 of table 808 and by storing the pre-calculated soft values: Assuming that the inphase and the quadrature components 
of the received signal each are quantised to n bits, the soft values of table 808 may be precalculated and tabulated for 
every different inphase and quadrature value, thereby further decreasing the required calculations, as the scaling and 
adding of step 805 are not necessary in this embodiment. However, such a table of pre-calculated soft values increases 
the memory consumption. Above, a=IRe(r)l and b=llm(r)l were defined as the absolute values of the real and imaginary 
25 parts of r, respectively, i.e. without sign information. Hence, a and b, each are represented by n-1 bits. If each of the 
pre-calculated soft values is to be represented by m bits, the total memory consumption of a full table is m.2 n " 416 
bits (for each of the 1 6 decision regions and each of the 4 bits, 2"-i different soft values are stored, each with a precision 
of m bits). For example, for n=4, the total memory is 512m bits. Note, however, that in this embodiment, the pre- 
calculated table still needs to be multiplied with the factor K. 
30 [0067] It is noted, that the above memory consumption may be further reduced by utilising the fact that many of the 
entries of table 808 are identical and by utilising the symmetry of the constellation of fig. 7. 

[0068] It is noted that the invention was described in connection with soft values defined as a log-likelihood ratio 
indicatingareliabilityvalueforthe bit values of areceived sequence. However, other definitions of soft values depending 
on the distance of the received signal to the signal points may be used as well. 
35 [0069] It is further noted that the signal constellations of figs. 2 and 7 are merely used as examples. The calculation 
of soft values according to the invention is not limited to these signal constellations. 

[0070] Finally, it is noted that the embodiments described in connection with figs. 3-6 are particularly well suited for 
large signal constellations, as they save memory, whereas the embodiment of figs. 8-9 is particularly well suited for 
implementing a medium-size signal constellation, e.g., 16QAM, as the save computational resources. 
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Claims 



1. A method of generating a reliability value (L k m ) for a received multilevel signal (r) in relation to a number of pre- 
*s determined signal symbols (S 1( .... S M ) each associated with a corresponding bit sequence including a first bit 

position (m); the reliability value being indicative of likelihood Information of receiving said multilevel signal; 
characterised in that the method comprises the steps of 

identifying (302; 502; 802) a first one ($ k ) of the number of signal symbols as being closest to the received 
so multilevel signal; and 

estimating (305, 306, 307; 506, 507; 804, 805) the reliability value based on stored information related to a 
second one of the number of signal symbols being closest to the first signal symbol and corresponding 
to a different binary value at the first bit position of the respective associated bit sequence than the first signal 
symbol. 



2. A method according to claim 1 , characterised in that the stored information comprises an identification of the 
second signal symbol, and the step of estimating the reliability value further comprises the steps of 
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determining (303) a first distance (5<,) between the received signal and the first signal symbol; and 
determining (305) a second distance (6 2 ) between the received signal and the second signal symbol. 

3. A method according to claim 1 , characterised in that the stored information comprises a second distance (5 3 ) 
5 between the first signal symbol and the second signal symbol, and the step of estimating the reliability value further 

comprises the step of determining (503) a first distance (5^ between the received signal and the first signal symbol. 

4. A method according to claim 2 or 3, characterised in that the step of estimating the reliability value comprises 
the step of determining (306, 307; 506, 507) a polynomial function of the first distance and the second distance, 

10 multiplied by a predetermined constant (K). 

5. A method according to claim 4, characterised in that the predetermined constant is selected depending on the 
noise distribution of the received multilevel signal. 

is 6. A method according to claim 1 , characterised in that the stored information is indicative of a number of functional 
relations between the received multilevel signal and the reliability value, and the step of estimating the reliability 
value further comprises the step of selecting 804 a functional relation of said number of functional relations de- 
pendant on the first signal symbol and the first bit position. 

20 7. A method according to claim 1 , characterised in that the stored information comprises for each signal symbol 
and bit position an approximation of the corresponding reliability value. 

8. A method according to any one of the claims 1 through 7, characterised in that the stored information is stored 
in a look-up table (308; 508; 808) indexed by the number of signal symbols and the bit positions of the number of 

25 bit sequences. 

9. A method according to any one of the claims 1 through 8, characterised in that the method further comprises.the 
step of providing the reliability value as an input to a decoder (106). 

30 10. A method according to any one of the claims 1 through 9, characterised in that the first signal symbol is identified 
by means of a slicer. 

11. A method according to any one of the claims 1 through 10, characterised in that the likelihood information com- 
prises a log-likelihood ratio. 

12. A method according to any one of the claims 1 through 11 , characterised in that the step of identifying the first 
signal symbol as being closest to the received multilevel signal comprises the step of identifying the first signal 
symbol as being closest to the received multilevel signal with respect to a Euclidean distance measure in a signal 
space. 

13. A method according to claim 12, characterised in that the signal space is related to the complex plane in quad- 
rature amplitude modulation. 

14. A method according to any one of the claims 1 through 13, characterised in that the number of signal symbols 
45 are associated with the number of bit sequences such that the bit sequences associated with all nearest neighbours 

of each signal symbol only differ from the bit sequence of that signal symbol at one bit position. 

15. An arrangement for generating a reliability value (L k m ) for a received multilevel signal (r) in relation to a number 
of predetermined signal symbols (S 1 ,...,S M ) each associated with a corresponding bit sequence including a first 

50 bit position (m); the reliability value being indicative of likelihood information of receiving said multilevel signal; 

characterised in that the arrangement comprises 

first processing means (104) adapted to identify a first one of the number of signal symbols as being closest 
to the received multilevel signal; 
55 - storage means (1 05) adapted to store information related to the first signal symbol and a second one of the 

number of signal symbols being closest to the first signal symbol and corresponding to a different binary value 
at the first bit position of the respective associated bit sequence than the first signal symbol; and 
second processing means (1 04) adapted to estimate the reliability value on the basis of the stored information. 
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16. An arrangement according to claim 15. characterised in that the storage means is 

cation of the second signal symbol, and the second processing means is further adapted to determine a first 
distance between the received signal and the first signal symbol and a second distance between the received 
signal and the second signal symbol. 

17. An arrangement according to claim 15. characterised in that the storage means is adapted to store a second 
distance between the first signal symbol and the second signal symbol, and the second processing means isfurther 
adapted to determine a first distance between the received signal and the first signal symbol. 

io 18 An arrangement according to claim 1 6 or 1 7. characterised in thatthe second processing means is further adapted 
to deterTne a polynomial function of the first distance and the second dfctance, multiplied by a predetermined 
constant. 

1 9. An arrangement according to claim 1 8, characterised in that the predetermined constant is selected depending 
15 on the noise distribution of the received multilevel signal. 

20. An arrangement according to claim 15, characterised in that the storage means is adapted to store information 
MtaM or a number of functional relations between the received multilevel signal and the reliability value, and 
the second processing means is further adapted to select a functional relation of said number of functional relations 

20 dependant on the first signal symbol and the first bit position. 

21. An arrangement according to claim 15, characterised in that the stored information comprises for each signal 
symbol and bit position an approximation of the corresponding reliability value. 

25 22. An arrangement according to any one of the claims 15 through 21 , characterised in that ^^^^ 
adapted to store the stored information as a look-up table indexed by the number of signal symbols and the bit 
positions of the number of bit sequences. 

23. An arrangement according to any one of the claims 15 through 22, characterised in that the first processing 
30 means further comprises a slicer. 

24. An arrangement according to any one of the claims 1 5 through 23, characterised in that the likelihood information 
comprises a log-likelihood ratio. 

35 25. An arrangement according to any one of the claims 15 through 24, characterised in that ^ P™^£ 
means is further adapted to identify the first signal symbol as being closest to the received multilevel signal with 
respect to a Euclidean distances in a signal space. 

26. An arrangement according to claim 25, characterised in that the signal space is related to the complex plane in 
40 quadrature amplitude modulation. 

27. An arrangement according to any one of the claims 15 through 26, characterised in that the number of signal 
symbols are associated with the number of bit sequences such that the bit sequences associated with all nearest 
neighbours of each signal symbol only differ from the bit sequence of that signal symbol at one bit position. 

28. A device (103) for receiving multilevel signals comprising an arrangement according to any one of the claims 15 
through 27. 

29. A device according to claim 28, characterised in that the device further comprises a decoder (106) adapted to 
so receive an input signal from the arrangement indicative of the determined reliability value. 

30. A device according to claim 28 or 29, characterised in that the device is a mobile terminal. 
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